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Abstract 

We investigate geometric phase (GP) of a qubit symmetrically coupled 
to a XY spin chain with three-spin interaction in a transverse magnetic field. 
An analytical expression for the GP is found in the weak coupling limit. It is 
shown that the GP displays a sharp peak or dip around the quantum phase 
transition (QPT) point of the spin chain. Without the three-spin interaction, 
the GP has a peak or dip around the critical point A = 1. If the three-spin 
interaction exists, the peak or dip position is obviously shifted away from 
the original position. This result reveals that the GP may be taken as an 
observable to detect both the existence and strength of multi-spin interaction 
in a spin chain. 
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1. Introduction 

Geometric phases (GPs) have been proposed as a typical mechanism for a 
quantum system to keep the memory of its evolution in Hilbert space. Since 
discovered first by Berry GPs have become an object of intense research 
in both the theoretical and experimental realm jij. The original concept of 
GPs has been generalized to the non-adiabatic case [jjj], the noncyclic case j3] 
and the degenerate state [^j. It is also noticed that GPs can be related to a 
number of important phenomena in physics 0] such as the Aharonov-Bohm 
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effect p]], the quantum Hall effect [8f and even applications in the quantum 
information processing 0,0, HQ- 

In recent years, relation between GPs and quantum phase transitions 
(QPTs) |13j in various closed many-body systems has attracted many inter- 



ests [14J, [15|, lis lHj • Carollo et al. [14( showed that GPs are much sensitive 
to controlling parameters of spin chains and can be exploited as a tool to 
detect critical regions of the systems. Zhu 15] investigated GPs in the XY 
spin chain and showed that the ground-state GP obeys the scaling behavior 
in the vicinity of QPTs. And it is also shown that the connection of GPs 
with the typical features of QPTs such as the scaling behavior, critical expo- 
nents and so on is not restricted to the XY spin chain model but universal 



for quantum many-body systems 15 



Since quantum systems are unavoidably to interact with their surround- 
ings, the time evolution of their states is generally nonunitary. Therefore, it 
is desirable to extend the concept of GPs from closed quantum systems to 
open ones. Until now, many approaches have been proposed for this purpose 



18, 19, 20, 21, 22 



The first work concerned with this point was given in 
Ref. 19], where GP is investigated purely as a mathematical problem. Based 
on the experimental context of quantum interferometry, Sjoqvist et.al 20 



introduced a definition of GP for mixed states undergoing unitary evolution. 



Tong and coworkers [22] developed a kinematic approach to GP for open 
quantum systems in nonunitary evolution led by the environments. Accord- 
ing to the definition of GP for open systems, many works concerned with the 
correction of environments to GPs have been done 23, 24, 25, 26, 27 . Yuan 



et.al |24| studied GP of a qubit coupled to an antiferromagnetic spin and 
found that the GP changes abruptly to zero when the spin chain undergoes 
a spin-flop transition. Lombardo et.al 25j computed the GP of a spin-1/2 
linearly coupled to a harmonic oscillator reservoir at arbitrary temperatures 
and estimated the time scale for experimentally measuring the GP. Recently, 
Villar et.al [26] investigated GP of a spin-1/2 particle in the presence of a 
composite environment, composed of an external bath and another spin-1/2 
particle. Their results show that the initial entanglement enhances the stur- 
diness of GP to decoherence. Cuchietti et.al 27J reported a measurement of 
GP for a spin-1/2 undergoing nonunitary evolution induced by coupling of 
an environment with a NMR quantum simulator. 

In previous studies of GPs with s pin environments, only is the nearest- 
neighbor spin interaction considered [24J, |28|. Besides the two-spin kind of 



interactions, however, multi-spin interactions may also exist in spin chains 



2 



29, 30, 31, 32 



As a result, phase transition points may be changed by 
the multi-spin interaction. In previous investigations, it has been shown 
that GPs of a qubit coupled to a spin chain change greatly around phase 
transition points of the spin chain and can signal out the appearance of 
phase transitions. According to this point, GPs are expected as a detector 



of phase transitions of quantum many-body systems [14j, [15|, [16|, [17| . With 
the same reason, we expect that GPs of a qubit coupled to a spin chain with 
multi-spin interaction may be used as an observable to detect the multi- 
spin interaction. In the present work, we shall investigate GPs of a qubit 
symmetrically coupled to an anisotropic XY spin chain, in which a three- 
spin interaction is included besides the nearest-neighbor spin interaction. 



According to the definition of GP given in Ref . [22] , we obtain an analytical 
expression for the GP of the central qubit in the weak coupling limit. Our 
results show that the variation of GP at the critical points of quantum phase 
transitions is much sensitive to the three-spin interaction, and the GP can be 
taken as a tool to detect the existence and strength of multi-spin interaction 
in spin chains. 

The present paper is organized as follows. In Section 2, the model is 
introduced. In Section 3, an analytical expression of the GP is obtained in 
the weak coupling limit and the effect of the multi-spin interaction on the 
GP is investigated. Finally, a brief summary is given in Section 4. 



2. The Model 

The model under consideration is composed of a central spin-1/2 (or 
qubit) and a iV-spin-1/2 chain. The central spin is symmetrically and trans- 
versely coupled to the circle spin chain, in which besides the nearest-neighbor 
spin interaction a three-spin interaction is also included. Meanwhile, a trans- 
verse magnetic field is homogeneously applied to each spin of the chain. The 
associated Hamiltonian reads 

N 

H = Wo-E 

i=i 

N N 

-a E (vti°t°i+i + °i-i°ttf+i) ~ K °i> W 
i=i i=i 

where rj is the transition frequency between the ground state (\g)) and the 
excited state (|e)) of the central qubit, 7 and a are the anisotropic parameter 



(1+7) 



-a, a 



(1-7) 



1 u i+i 



3 



of the two-spin interaction and the strength of the three-spin interaction, 
respectively, A is the coupling constant of the spin chain with the transverse 
magnetic field, and g is the coupling strength between the central qubit and 

x, y, z) are the Pauli matrices for spin at the 
e) (e| — \g) (g\) is the population inversion 



the spin circle. In ([!]), <7[™(m 



Ith site of the spin chain and <Tq( 
operator for the central spin. 

The free-motion term of the central qubit in ([T]) can be removed by the 
rotating transformation exp (— iria^i). Then, following the approach pro- 



posed in Ref. [13l . |32|, the rotated Hamiltonian can be fully diagonalized. 
Since [ctq,^" 1 ] = 0, an operator-valued parameter A = A + ga^ is a con- 
served quantity. Thus A can be treated number during diagonalizing 
the Hamiltonian. Obviously, A has two eigenvalues A n = A + (— l) n g with 

71 = 0,1. 



By introducing the Jordan-Wigner transformation 13 



en 



n 

i<l 



1 - 2ck 



C'l 



af 



i<l 



- lc\ci) (ci-cj, . 



(2) 



(3) 



of = 1 - 2c\c h (4) 

where c\ and q are the mapped spinless and fermionic creation and annihi- 
lation operators. Substituting Eqs. (J2J)- (J4] ) into ([[]), one obtains 



A' 



H 



i=i 

+2a (c\_ lCl+1 + 4 +1 c,_ x ) + A (l - 2cJq) } . 



(5) 



Then by means of the Fourier transformation [13 



k=-M 
M 



(6) 



(7) 



k=-M 
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with M = N/2 for N even and M = (N - 1) /2 for N odd, one gets 
if = 2 (A — cos A;a — 2a cos 2A;a) d].c4 + 27 sin fca ^d£<it fc — d- k dk^ , (8) 



£:>0 



where d k and dj, are the fermionic annihilation and creation operators in the 
momentum space, respectively. 

In order to diagonalize the Hamiltonian (jSJ), we introduce the Bogoliubov 



transformation 13 



qA qA 

7fe,A = cos -y d fc - i sin y d! fe . (9) 

By substituting ([2]) into (jSj), the Hamiltonian ([T]) can be fully diagonalized 
and written in the form 



H = J2^k,A (7l,A7fc,A-^J 



where 



and 



(10) 



£lk,A = 2y (A — cos ka — 2a cos 2ka) + 7 2 sin ka, (11) 



, , 7 sin ka , 

9t = arctan j- — (12) 

A — cos ka — 2a cos 2ka 



with a = 2-n/N. 

It is noted that the fully diagonalized Hamiltonian for the pure spin chain 
can be obtained by setting g = in Eq. (ITUj) . Correspondingly, the energy 
spectrum Q kj \, the parameter 9% and the mode operator 7 fc A for the pure 
spin chain can be obtained just by changing A into A in Eqs. fl9|)-f lT2|) . It can 
be easily proved that both the mode operators 7^ for the spin chain and 
7k ( A for the qubit-chain coupled system are related by 

7fc,A = 7fe,A cos a kA - ryl fc A sin a kA (13) 

with a kA = (e£ - el) 12. 
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3. Geometric Phase of the Central Qubit 

Let us suppose that at t = the central qubit is decoupled from the spin 
environment, and the qubit is in a state 

|0(O)) = cos^|e) + sin^>. (14) 

The ground state of the spin chain is defined as 7^ \G) X = 0. It can be 

written in the form | G) x = J~[ (cos -|- \0) k |0}_ fc + i sin \l) k ), where 

fe>o ^ " ' 
|0) fe and |l) fc are the vacuum and single excitation states of the /cth pure en- 
vironment mode, respectively. Then according to (|T3|) . \G) X can be rewritten 

as 

\°)x = II ( C0S ak > A ~ i Sil1 ^^Ia^-^a) I G )a > ( 15 ) 
fc>0 

where |G) A is the ground state of the Hamiltonian (fTOl) and satisfies 7^ \ G) A = 
0. 

At time t, the qubit-environment coupled system unitarily evolutes into 
the state p(t) = U (t) (\<p (0)) ® |G) AA (G| ® (0 (0)|) f/ 1 " (t) with the time evo- 
lution operator U (t) = exp (—iHt). By tracing p (t) over variables of the 
spin chain, we obtain the reduced density matrix for the central qubit 



cos 2 I 



±sin/3F(t) 



* ( ' )= U^W) J' (16) 

where F (t) is the decoherence factor and is given by 

F(t)= x (G\Ul(t)U Ao (t)\G) x . (17) 

The evolution operators U\ (t) and U Al (t) are obtained from U (t) by re- 
placing A with X + g and A— g, respectively. From Eqs. (fT5| and ( ITT)) , we can 
work out the explicit expression for the modulus of the decoherence factor 



= ]J{AB cos (Q k>Ao -n k>Al )t-AB 2 sin 2 (a k , Ao -a k>Al ) (18) 

fc>0 

1 /2 

+ 1 - sin 2 (2o fciAo ) sin 2 {& k ,A Q t) - sin 2 (2a fe)Al ) sin 2 {Q k)Al t)} 
where A = sin (2a fc;Ao ) sin (2^^), B = 2 sin (O fe) A t) sin (f^A^)- 
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According to the GP definition for an open system, which is proposed in 
Ref. [22J, the GP acquired by the qubit in a quasi period is given by 



$ = arg V y/e k (0) e k (T) (e k (0)| e k (T)> c -J?M<)|fl/*k*«>>* j (19) 



where e k (t) and \e k (£)) are the eigenvalues and the corresponding eigenstates 
of the reduced density matrix ffTB]) . and T = 2ir/r] is the time evolution cycle 
of the qubit when it is isolated from the environment. If the qubit is coupled 
to the spin chain, the evolution of its state is no more periodic. However, if 
the coupling is weak, we may consider a quasi cyclic path P : t G [0, T]. 

The instantaneous eigenvalues and the corresponding eigenvectors of ([16]) 
are found to be 

e± (t) = ~ (l ± ^cos^ + sin 2 /?!^)! 2 ) (20) 
|e ± (t)> = e"^ sin Ml |e> + cos M> \g) , (21) 

with 

cot ± J cot 2 /3+ \F(t)\ 2 
(3 ± (t) = 2 arctan \F(tj\ ' (22) 

Since e_ (0) = 0, Eq. fTl9l) shows that only the eigenvalue e + (£) and corre- 
sponding eigenvector \e + (£)) have contribution to the GP of the qubit. Upon 
substituting Eqs. ([IHD-flZSD into Eq. (TBI]), we obtain the GP of the central 
qubit 

j>2n/ri a 

§ = r] dt sin 2 — . (23) 



o 



2 

If the qubit has no interaction with the spin chain, i.e. g = 0, the well- 
known result $ = 7i (1 + cos 0) can be recovered from (|23|) . When the qubit 
is coupled to the spin chain, the GP will be modified. In order to get a basic 
and clear intuition about the effect of the three-spin interaction on the GP 
acquired by the qubit in a quasi cycle, we first make an approximate analysis 
on the GP. 

From Eqs. ( 122]) and f l23|) . we see that the correction of the spin chain to 
the GP is completely included in the decoherence factor F (t). Therefore, we 



first analyze the decoherence factor. For this purpose, following Refs. [33 
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we introduce a cutoff number K C} which determines the largest energy scale 
of the spin chain, and define the partial product 



K c 



\F(t)\ c = l[F k >\F(t)\, (24) 



k>0 



and the corresponding partial sum S (t) = \n\F (t)\ c = ^2 k ° Q \hiFk\. For 
small k, we can expand (fTTj) as a power series of k and obtain 

n kAn = 2 |A n - 1 - 2a\ (25) 

and 

(-l) n+1 27T 7 % , , 

Smat ^ JV|(A„-l-2 a )(A-l-2a)r (26) 
• / \ —2n^kg 

Sill [Oth A n — Oth Ai ) ~ r~r~: —7- rr. (27) 

v *° ' lJ iV|(A -l-2a)(Ai-l-2a)| 
In this way, the partial sum S (t) can approximately be written as 

^LVP. (28) 

2iV2(A-l-2a) 2 ^ 

If fT c is small, in the thermal dynamical limit iV — > oo, one may ignore all 
the terms related to N~ l with I > 3 in (1281) . Consequently, in a short time, 
we have 

\F(t)\ c ^e- Tt \ (29) 

where r = 8E (K c ) 7 V / (A - 1 - 2a) 2 and E (K c ) = 4ir 2 K c (K c + 1) x 
(2K C + 1) / (QN 2 ). This result shows that the decoherence factor would de- 
cay in a Guassian-type behaviour in a short time. Compared with that in 
the Ising spin environment, the decay rate r is modulated by the three-spin 



interaction a of the chain 33 



Upon substituting f )29|) into fl22|) . keeping all the terms up to the second 
order of the coupling strength g and completing the integral ( 12"3"|) . we obtain 
the approximate expression for the GP 

^ 64£(K c )7r 3 7 2 oos/3sin 2 /3 

« = + + VA-l-*.)' 9 <30) 
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In the above expression, the first term comes out when the qubit is not 
coupled to the spin chain and undergoes a unitary evolution. The second 
term is the modification induced by the spin chain. Since the modification is 
proportional to the squared coupling strength g, and the squared reciprocal 
of the transition frequency rj and the criticality factor (A — 1 — 2a) ~ 2 , one 
may expect that the three-spin interaction can dramatically change the GP 
modification around the point A — 1 — 2a = except ft = 0,7r/2 and ir. 
Meanwhile the GP modification is positive when ft < 7r/2 and negative when 
ft > 7r/2, and changes sign around ft = ir/2. 

In order to get the entire picture about the effect of the three-spin inter- 
action on the GP, we have numerically investigated the variation of GP with 
various parameters of the spin chain according to Eq. (I2"5j) . 




Figure 1: The GP versus the three-spin interaction strength a for the different initial states 
in the Ising limit (7=1) and without the external field (A = 0). The other parameters 
are chosen as N = 501, g = 0.03 and 77 = 2tt/3. 

When the transverse magnetic filed is absent, i.e. A = 0, Eq. ( 1301) shows 
that the GP modification may dip or peak around the point a = —0.5, 
depending on the initial state of the qubit. In Fig. 1, the GP is plotted as a 
function of the three-spin interaction strength a for the different initial states 
in the Ising limit (7 = 1). The features of the GP on the parameter ft shown 
in Fig. 1 are same as expected from Eq. ( l30l) . For a fixed value of ft, except 
the points ft = 0,7r/2 and n, the GP displays a peak or dip at the critical 



9 



point a = —0.5. Moreover, it is noted that the GP also has a peak or dip at 
the critical point a = 0.5 which can not be expected from Eq. (130 p . In Eq. 
(18), the variable k takes values from to N/2. Eqs. (11)-(12) depend on the 
variable k through the trigonometric functions sm(2nk/N) and cos(27r/c/iV). 
Thus, Eqs. (11)-(12) must have dips or peaks around k = N/2 if it has dips 
or peaks around k — 0. It is this reason that leads to the appearance of 
peaks or dips of the GP at a = 0.5. Fig. 1 also shows that the critical point 
position does not change with variation of /3, and is determined only by the 
three-spin interaction. 




Figure 2: The GP versus the three-spin interaction a with (a) different vaiues of of the 
transition frequency r\ and N — 501, g — 0.03, j3 = tt/5 and 7 = 1; (b) different values of 
the qubit coupling strength g and N = 501, 77 = n/5, (3 = it/5 and 7 = 0.4; (c) different 
values of the spin chain size N, and g = 0.03, r\ — 2tt/3, (3 = 7r/3, 7 = 1 and (d) different 
values of the anisotropic parameters 7, and g — 0.03, N — 501, rj = tt/5, (3 — n/5. 

In Figs. 2, the GP is plotted as a function of a with various values of 
the transition frequency 77, the coupling strength g and the spin chain size 
N for the case of A = 0, respectively. The GP versus a is shown in Fig. 
2(a). It is observed that the variation of the GP is sharpened around a = 0.5 
by decreasing the transition frequency rj, which is in accordance with the 
analytical result of Eq. ( 1301) . The influence of the coupling strength g on 
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the GP is shown in Fig. 2(b). Obviously, the peaks of the GP become less 
pronounced with increasing g. It means that as a tool of detecting QPTs 
the GP is powerful only in the weak coupling case. The influence of the spin 
chain size N on the GP is given in Fig. 2(c). It is seen that the peak becomes 
less sharp with increasing N. The GP of the qubit coupled to the XY spin 
chain(7 ^ 1) is shown in Fig. 2(d). Obviously, the larger the anisotropy 
parameter is, the more sharp the GP peak is. 

Combing Fig. 1 with Fig. 2, we see that by varying the parameters the 
peak values of the GP is changed but the peak position of the GP is retained. 
Thus, the QPT at |a| = 0.5 induced by the intrinsic three-spin interaction 
cab be clearly singled out by the GP of the central qubit. 




Figure 3: The GP as a function of A with different values of the three- spin interaction 
strength a in the Ising limit. The other parameters are chosen as7 — 1, N = 501, g = 0.03, 
T) = 2n/3 and (i = tt/3. 

In Fig. 3, the GP is shown as a function of A with various values of 
the three-spin interaction strength a in the Ising limit. When the three-spin 
interaction is absent, i.e. a = 0, the spin chain has a ground-state phase 
transition at A = 1. As shown in Fig. 3, the GP displays a peak at A = 1 
when a = 0. When the three-spin interaction is switched on (a ^ 0), as 
expected from ( 1301) . the three-spin interaction changes the original critical 
point and the GP peak has a 2a displacement from the point A = 1. 

In Figs. 4, the GP is plotted as a function of A with a = 0.2 for different 
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Figure 4: The GP versus the coupling constant A for a = 0.2 with (a) different values of of 
the transition frequency 77 and N = 501, g = 0.03, f3 = it/5 and 7 = 1; (b) different values 
of (3 and N = 501, 77 = 2tt/3, g = 0.03 and 7 = 1; (c) different values of the spin chain 
size N, g — 0.03, r\ = 2ir/3, f3 = n/3, 7 = 1 and (d) different values of the anisotropic 
parameters 7, and g = 0.03, N — 501, ij = it/5, f3 = tt/5. 
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values of the transition frequency 77, the initial state coefficient /3, the spin 
chain size N and the anisotropy parameter 7. In the Ising limit, by setting 
N = 501, g = 0.03, 7 = 1 and /3 = 7r/5, the variation of the GP is shown 
against A in Fig. 4(a). It is observed that the GP peaks have a 0.4 displace- 
ment from A = 1 and become more sharp around A = 1.4 with decreasing 
77, much similar to that in Fig. 2(a). By choosing N = 501, 77 = 2n/3 and 
g = 0.03, Fig. 4(b) shows the influence of /3 on the GP. It can be seen that 
the GP peaks also move to the point A = 1.4 but the peak value decreases 
with increasing /3. The influence of the spin chain size N on the GP is shown 
in Fig. 4(c) with the parameters g = 0.03, 77 = 2n/3 and (5 = ir/3. It is seen 
that the GP peak value increases with increasing N. In Fig. 4(d), the GP is 
plotted as a function of the anisotropy parameter with N = 501, 77 = 7r/5 and 
(3 = 7r/5. It is clear that the GP peak becomes pronounced but the critical 
point is not shifted as the anisotropy parameter increases. 

4. Summary 

We investigate geometric phase (GP) of a qubit symmetrically coupled 
to an anisotropic XY spin chain with three-spin interaction in a magnetic 
field. An analytical expression for the GP is obtained in the weak coupling 
limit. We find that the GP displays a peak or dip at the quantum phase 
transition (QPT) points of the spin chain, depending on the initial state 
of the qubit. Without the three-spin interaction, the QPT appears at the 
critical point A = 1 and the GP has a peak or dip around this critical point. 
If the three-spin interaction exists, the peak or dip position of the GP is 
obviously shifted away from the original position (A = 1). Moreover, it is 
analytically and numerically shown that the position shift of the GP peak or 
dip is determined only by the three-spin interaction. Thus, it comes to the 
conclusion that the GP can be used as a tool to detect both the existence 
and strength of multi-spin interaction in a spin chain. 
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